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Abstract. We present a new set of accurate formulae for the computation of random errors in the measurement of 
atomic and molecular line-strength indices. The new expressions are in excellent agreement with numerical simulations. 
We have found that, in some cases, the use of approximated equations can give misleading line-strength index errors. It 
is important to note that accurate errors can only be achieved after a full control of the error propagation throughout 
the data reduction with a parallel processing of data and error images. Finally, simple recipes for the estimation of 
the required signal-to-noise ratio to achieve a fixed index error are presented. 
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1. Introduction 

Line-strength indices have proven to be excellent tools for 
the study of absorption features in the spectra of astro¬ 
nomical objects. In particular, since the pioneering work of 
Faber (1973), the behaviour of indices in composite stellar 
systems has supplied fundamental clues for understand¬ 
ing the history of their stellar populations (e.g. Burstein 
et al. 1984; Bica et al. 1990; Gorgas et al. 1990; Worthey 
et al. 1992; Bender et al. 1993; Jones & Worthey 1995; 
Davies 1996). At the light of up-to-date stellar popula¬ 
tion models, line-strength indices can provide constrains 
about important parameters such as mean age, metallic- 
ity, or abundance ratios (e.g. Worthey 1994; Vazdekis et 
al. 1996; Bressan et al. 1996). In addition, line-strength in¬ 
dices have been widely employed to determine basic stellar 
atmospheric parameters (Rich 1988; Zhou 1991; Ternclrup 
et al. 1995) and star cluster abundances (Brodie & Hanes 
1986; Mould et al. 1990; Gregg 1994; Minniti 1995; Huchra 
et al. 1996), as well as to confront cosmological issues by 
providing new distance indicators (Dressier et al. 1987) 
and a powerful tool to investigate the redshift evolution 
of galaxies (Hamilton 1985; Chariot & Silk 1994; Bender 
et al. 1996). 

Obviously, the suitability of line-strength indices to in¬ 
vestigate the above items relies on a proper determination 
of the associated index errors. For illustration, it is inter¬ 
esting to note how errors in the measurement of absorption 
features translate into uncertainties in the derived mean 
age and metallicity of an old stellar population. Taking the 


predictions of the single-burst stellar population models of 
Worthey (1994) as a reference, typical errors of 0.20 A in 
the Lick indices Fe5270 and H (3 (see below) translate into 
A[Fe/H] ~ 0.18 dex and Aage ~ 6 Gyr respectively (for 
a composite population of 10 Gyr and metallicity around 
[Fe/H]=0.3). Another example in which it is extremely 
important to obtain reliable index errors is the analysis of 
the intrinsic scatter of relations such as that of Mg 2 with 
velocity dispersion in elliptical galaxies (Schweizer et al. 
1990, Bender et al. 1993). 

An accurate error determination is therefore needed 
in order to draw confident interpretations from observed 
data. Although this requires the estimation of both ran¬ 
dom and systematic errors, in this paper we will concen¬ 
trate in the former. Whereas random errors can be read¬ 
ily derived by applying statistical methods, systematic er¬ 
rors do not always allow such a straight approach. Sev¬ 
eral authors have dealt with the most common sources 
of systematic errors in the measurement of line-strength 
indices. These include: flux calibration effects (Gonzalez 
1993; Davies et al. 1993; Cardiel et al. 1995), spectral res¬ 
olution and velocity dispersion corrections (Gorgas et al. 
1990; Gonzalez 1993; Carollo et al. 1993; Davies et al. 
1993; Fisher et al. 1996; Vazdekis et al. 1997), sky sub¬ 
traction uncertainties (Saglia et al. 1993; Davies et al. 
1993; Cardiel et al. 1995; Fisher et al. 1995), scatter light 
effects (Gonzalez 1993), wavelength calibration and ra¬ 
dial velocity errors (Cardiel et al. 1995; Vazdekis et al. 
1997), seeing and focus corrections (Thomsen & Baum 
1987; Gonzalez 1993; Fisher et al. 1995), deviations from 
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linearity response of the detectors (Gorgas et al. 1990; 
Cardiel et al. 1995), and contamination by nebular emis¬ 
sion lines (Gonzalez 1993; Goudfrooij & Emsellem 1996). 

Although it is always possible to obtain estimates of 
index errors by performing multiple observations of the 
same object, analytical formulae to avoid such an observ¬ 
ing time effort are clearly needed. In this paper we present 
a set of analytical expressions to derive reliable errors in 
the measurement of line-strength indices (hereafter we use 
error to quote random errors exclusively). General index 
definitions are given in Section 2. Section 3 describes how 
a proper treatment of error propagation throughout data 
reduction is a prerequisite to compute confident index er¬ 
rors. Random errors can be obtained through numerical 
simulations, presented in Section 4, or analytically. Previ¬ 
ous works are briefly reviewed in Section 5, with special 
emphasis to the approach followed by Gonzalez (1993). 
Our final set of formulae is presented in Section 6. Sec¬ 
tion 7 gives some recipes to estimate the required signal- 
to-noise ratios to achieve a fixed index error. 

2. Definitions 


of broad bandpasses with the continuum regions located 
far from the central feature. On the other hand, atomic 
indices, which measure the absorption of atomic spectral 
features, have narrower and neighboring bandpasses. 

Throughout this paper we use the definitions given by 
Gonzalez (1993), in which atomic (/ a ) and molecular (J m ) 
indices are defined as follows: 


/a = 


(1 — 5(A)/C(A)) dA 


/ m = —2.5 log 


A7 S(X)/C(X) dX 


10 ' 


Ac 2 A Cl 


( 3 ) 

( 4 ) 


where A Cl and A C2 are the limits of the central bandpass 
(in A). The local pseudo-continuum C(A) is derived by 


C(A) = Sb+ 5 r --where (5) 

/X'T” Aj) /Af /A£) 


s b = Xbl 


ft S(X) dA 


(A b 2 Ae,i ) 


S r = 


f£ s W jA 

(Xr 2 A ri ) 


( 6 ) 


The strengths of absorption spectral features have been 
measured in different ways so far. However, although with 
slight differences among them, most authors have em¬ 
ployed line-strength indices with definitions close to the 
classical expression for an equivalent width: 

Wa(A) = f (1 - 5(A)/C(A)) dA (1) 

J line 

where 5(A) is the observed spectrum and C(A) is the lo¬ 
cal continuum usually obtained by interpolation of 5(A) 
between two adjacent spectral regions (e.g. Faber 1973; 
Faber et al. 1977; Whitford & Rich 1983; Gorgas 1987; 
Brodie & Huchra 1990; Gonzalez 1993; Rose 1994). In or¬ 
der to avoid subjective determinations of the local con¬ 
tinuum, line-strength indices following Eq. ([]]), referred 
as atomic indices , are completely characterized by three 
wavelength regions (bandpasses). The spectral feature of 
interest is covered by the central bandpass whereas the 
other two bandpasses, located towards the red and blue of 
the central region, are employed to define the continuum 
reference level through a linear interpolation. As pointed 
out by Geisler (1984) (see also Rich 1988), at low spectral 
resolution a pseudo-continuum is measured instead of a 
true continuum. 

Line-strength indices are sometimes measured in mag¬ 
nitudes using: 

/(mag) = -2.5 log 10 ^1 - (2) 

where AA is the width of the central bandpass. These 
line-strengths, referred as molecular indices since they are 
used for molecular-band features, are defined with the help 


A b = (Ab 1 + Ab 2 )/2, 


X r = (A ri + X r2 )/2 (7) 


being X bl , A& 2 , A ri , and X r2 the limits of the blue and red 
bandpasses respectively. 

Although simplified versions of these expressions (i.e. 
considering a constant continuum flux along the central 
bandpass, or replacing the integrals by mean values) yield 
similar results at intermediate resolution, the more ac¬ 
curate formulae must be favoured in order to guarantee 
the comparisons with high-resolution high-S/N (signal-to- 
noise ratio) spectra (specially for asymmetric indices, as 
already noted by Worthey et al. 1994). 

Probably, the most widely index definition system 
employed so far is that established by the Lick group 
(Burstein et al. 1984, 1986; Faber et al. 1985; Gorgas et al. 
1993; Worthey et al. 1994). In Table 1 we list the exact def¬ 
initions (as given by Trager 1997) for the 21 indices which 
constitute the extended Lick system (see also Gonzalez 
1993). This table also includes the definitions given by 
Diaz et al. (1989) for the Call triplet in the near-infrared. 
In this paper we will concentrate on the analysis of line- 
strength errors for these particular indices, although the 
derived analytical expressions are valid for any general in¬ 
dex following Eq. (||) or (||). 

Another interesting spectral feature which will be stud¬ 
ied in this paper is the amplitude of the A4000-A break 
(D 4 ooo)- We adopt here the definition given by Bruzual 
(1983): 


D4000 = 


r 4250 
J 4050 
r3950 
J 3750 


5(1/) dA 
S(v) dA 


( 8 ) 


This index can be considered as a pseudo-color, being the 
combination of 5(i/) and dX due to historical reasons. 
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Table 1 . Bandpass definitions for the Lick index system —using the revised bandpass limits as given by Trager (1997)—, 
together with the Call triplet according to Diaz et al. (1989). The S/N constants ci, C 2 and C 3 are explained in Section 7. 


Index Name 

Central Bandpass (A) 

Continuum Bandpasses (A) 

Ci 


Atomic Indices 


Cl 

C2 

Ca4227 

4222.250-4234.750 


4211.000-4219.750 

4241.000-4251.000 

4.604 

0.3684 

G4300 

4281.375-4316.375 


4266.375-4282.625 

4318.875-4335.125 

8.537 0.2439 

Fe4383 

4369.125-4420.375 


4359.125-4370.375 

4442.875-4455.375 

13.220 

0.2580 

Ca4455 

4452.125-4474.625 


4445.875-4454.625 

4477.125-4492.125 

7.038 

0.3128 

Fe4531 

4514.250-4559.250 


4504.250-4514.250 

4560.500-4579.250 

11.299 

0.2511 

Fe4668 

4634.000-4720.250 


4611.500-4630.250 

4742.750-4756.500 

17.757 

0.2059 

H/3 

4847.875-4876.625 


4827.875-4847.875 

4876.625-4891.625 

7.301 

0.2539 

Fe5015 

4977.750-5054.000 


4946.500-4977.750 

5054.000-5065.250 

16.455 

0.2158 

Mgb 

5160.125-5192.625 


5142.625-5161.375 

5191.375-5206.375 

8.032 0.2472 

Fe5270 

5245.650-5285.650 


5233.150-5248.150 

5285.650-5318.150 

9.250 

0.2313 

Fe5335 

5312.125-5352.125 


5304.625-5315.875 

5353.375-5363.375 

10.741 

0.2685 

Fe5406 

5387.500-5415.000 


5376.250-5387.500 

5415.000-5425.000 

7.256 

0.2893 

Fe5709 

5696.625-5720.375 


5672.875- 5696.625 

5722.875- 5736.625 

6.362 

0.2679 

Fe5782 

5776.625-5796.625 


5765.375-5775.375 

5797.875-5811.625 

6.134 

0.3067 

NaD 

5876.875-5909.375 


5860.625-5875.625 

5922.125-5948.125 

8.113 0.2496 

Cal 

8483.000-8513.000 


8447.500- 8462.500 

8842.500- 8857.500 

8.852 0.2951 

Ca2 

8527.000-8557.000 


8447.500- 8462.500 

8842.500- 8857.500 

8.330 0.2777 

Ca3 

8647.000-8677.000 

Molecular Indices 

8447.500- 8462.500 

8842.500- 8857.500 

7.750 0.2583 

C 3 

CNi 

4142.125-4177.125 


4080.125- 4117.625 

4244.125- 4284.125 

0.2241 

cn 2 

4142.125-4177.125 


4083.875-4096.375 

4244.125-4284.125 

0.2691 

Mgi 

5069.125-5134.125 


4895.125- 4957.625 

5301.125- 5366.125 

0.1662 

Mg 2 

5154.125-5196.625 


4895.125- 4957.625 

5301.125- 5366.125 

0.1933 

TiOi 

5936.625-5994.125 


5816.625- 5849.125 

6038.625- 6103.625 

0.1824 

Ti0 2 

6189.625-6272.125 


6066.625- 6141.625 

6372.625- 6415.125 

0.1568 
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3. Propagation of errors throughout data reduc¬ 
tion 

The aim of the reduction process is to minimize the influ¬ 
ence of data acquisition imperfections on the estimation 
of the desired astronomical quantity (see Gilliland 1992 
for a short review on noise sources and reduction pro¬ 
cesses of CCD data). For this purpose, one must perform 
appropriate manipulations with the data and calibration 
frames. The arithmetic work involved in this process must 
be taken into account in order to get reliable estimates of 
line-strength errors. 

In order to trace in full detail the error propagation, 
error frames must be created at the beginning of the re¬ 
duction procedure. After this point, error and data frames 
should be processed in parallel, translating the basic arith¬ 
metic manipulations performed over the data images into 
the error frames by following the law of combination of 
errors. 

The starting point is the creation of initial error images 
with the expected RMS variances from photon statistics 
and read-out noise. For a single spectrum: 

Aj} = - g N c [j}+a 2 RN [j} (9) 

where <J 2 [j] is the variance in the pixel [j] ( a\j ] in number 
of counts, ADU, —analogic to digital number—), g the 
gain of the A/D converter (in e _ /ADU), N c [j] the number 
of counts in the pixel [j] (after the bias- level subtraction), 
and ct R n[/] is the read-out noise (in ADU)El. 

Some of the reduction steps that may have a non neg¬ 
ligible effect in the index errors are flatfielding, geometri¬ 
cal distortion corrections, wavelength calibration, sky sub¬ 
traction and rebinning of the spectra. Note that if error 
spectra were computed from the final number of counts 
in the reduced data frame, index errors would tend to be 
underestimated. The extra benefit of a parallel processing 
of data and error frames is the possibility of obtaining, at 
any time of the reduction, the variation in the mean S/N 
ratio produced by a particular reduction step. Under these 
conditions, it is possible to determine which parts of the 
reduction process are more sensitive to errors and even 
to decide whether some manipulations of the data images 
can be avoided, either because the resulting S/N ratio is 
seriously reduced and the benefit of the product insignif¬ 
icant, or because such manipulations become a waste of 
time. 

A full error propagation in the reduction of spectro¬ 
scopic data has been previously implemented by Gonzalez 

1 Note that the apparent dimensional inconsistency of Eq. [] is 
not real, and arises from the fact that one of the properties 
of the Poisson distribution is that its variance is numerically 
equal to the mean expected number of events. 


(1993), and it is also included in_the reduction package 
ReD^E (Cardiel & Gorgas 1997)1. 

To obtain errors on line-strength indices from the re¬ 
duced data and error spectra, two different approaches 
can be followed. On one hand, analytical formulae to eval¬ 
uate index errors as a function of the data and variance 
values in each pixel can be applied (Gonzalez 1993). An¬ 
other method is to simulate numerically the effect of the 
computed pixel variances in the index measurements (e.g. 
Cardiel et al. 1995). Both techniques are examined in the 
next sections. 

4. Numerical simulations 

The effect of random noise in the spectra of astronomi¬ 
cal objects can be simulated by introducing in each pixel 
Gaussian noise computed as: 

7 Z[j] = x a[j] yj— ln(l — n) cos(27rr 2 ) (10) 

where <j 2 [j] is the variance in the pixel [j], and rq and 
r 2 two random numbers in the range rq,r 2 £ [0,1). Af¬ 
ter the creation of Af s j mu i synthetic spectra, index errors 
can be derived as the unbiased standard deviation of the 
Asimui measurements of each index (typically we employed 
N s imul ~ 1000). 

For illustration, we show the results of numerical simu¬ 
lations using a high S/N ratio spectrum of the bright star 
HR 3428 as a template. We have assumed this spectrum 
to be noiseless, i.e. the nominal values of the measured 
line-strength indices in this spectrum were considered as 
error free reference values. By dividing this template spec¬ 
trum by different constant factors, we built a set of syn¬ 
thetic error spectra which were used to derive index er¬ 
rors as a function of the mean (S/N-ratio)/A for different 
atomic and molecular indices (Fig. |l|) . In a logarithmic 
scale, there is a clear linear correlation between the esti¬ 
mated relative error and the (S/N-ratio)/A. In addition, 
and as it should be expected, at a fixed (S/N-ratio)/A rel¬ 
ative errors for atomic indices (with narrow bandpasses) 
are larger than for molecular indices. 

5. Analytical formulae: previous works 

Several authors have employed analytical formulae to es¬ 
timate errors in the measurement of line-strength indices 
(Rich 1988; Brodie & Huchra 1990; Carollo et al. 1993; 
Gonzalez 1993). However, the derived error expressions 
have always been obtained from approximated versions 
of Eqs. (|[) and (|j), in which integrals are replaced by 
averaged number of counts in the bandpasses. Although 
this approximation works properly under certain circum¬ 
stances, some additional assumed simplifications seriously 
constrain the suitability of such formulae. The expressions 

2 Available at: _ 

http: / / www.ucm.es/OTROS/Astrof/reduceme/reduceme.htrnl 
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Fig. 1. Relative errors from numerical simulations in the mea¬ 
surement of 21 line-strength indices in the bright star HR 3428, 
as a function of the mean S/N-ratio per A. 


= S 2 J J X " ^ dX (18) 

and d 2 (A) is the variance of S( A) at the wavelength A. 

In order to check the accuracy of Eqs. © and (|l2|) we 
compared the predictions of such expressions with the re¬ 
sults obtained from numerical simulations. We took spec¬ 
tra and associated error spectra of a large homogeneous 
sample of 350 standard stars from the Lick library (ob¬ 
served with RBS —Richardson Brealy Spectrograph— at 
the JKT Jacobus Kapteyn Telescope— of the Roque de 
los Muchachos Observatory, La Palma, February 1995). 
The error spectra were obtained by following a parallel 
reduction of data and error frames, as it has been previ¬ 
ously described. As it is apparent from Fig. there is an 
excellent agreement between Gonzalez’s formulae and the 
simulations. 


derived by Rich (1988) do not take into account readout 
noise (which can be important at low count levels) nor the 
effect of the sky subtraction. In the formulae employed 
by Brodie & Huchra (1990) and Carollo et al. (1993), al¬ 
though these factors are considered, important reduction 
steps such as flatfielding (in particular slit illumination) 
are not taken into account. All these problems are readily 
settled with a parallel reduction of data and error frames 
(Gonzalez 1993, section 3 this paper). Although this ap¬ 
proach requires a more elaborate reduction, it yields the 
most confident results since it allows the application of 
analytical formulae which consider the final error in each 
pixel of the reduced spectra. 

The most accurate set of analytic formulae published 
so far for the computation of errors in line-strength indices 
are those presented by Gonzalez (1993). As a reference, we 
reproduce here his equations: 

<*%] = 


Sc 

Cc 


where 


OSc 

S c 


Uk / A r - A e 
C 2 \ Ar - A(, 


A b A c 
Cl V A r - A b 


_ 2.5 x 10 0 - 4J " 

a[ini 2.3026(A C2 — A Cl ) CT a 


A c = (A Cl + A C2 )/2 
C c = C{ A c ) 


Sc = 


0 vsJSc ) 2 = V 


<4 = Si/ 


5(A) dX 
S 2 ( A) 

s 2 {x) 

°' 2 (A) 


dX 




dX 


' Abj. 


( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 
(17) 



Relative Error (simulations) 


Fig. 2. Comparison of the relative errors employing Gonzalez’s 
formulae with those obtained from numerical simulations. 
The symbols corresponds to the measurement of six different 
line-strength indices in a sample of 350 stars from the Lick 
library. 


Unfortunately, the showed concordance between the 
analytical formulae and the simulations can not be ex¬ 
trapolated to any general situation. In particular, one of 
the approximations used to derive Eqs. & and ( |l2| ) is 
the assumption that the mean signal-to-noise ratio in each 
bandpass can be computed as the quadratic sum of the in¬ 
dividual signal-to-noise ratio in each pixel (Eqs. |Tg 18). In 
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other words 


Eil p i xels S(X t) 
E t p f eIS ^ 2 (A.) 


pixels 

E 


i= 1 


s 2 {x l ) 

O ' 2 (A,) 


(19) 


where A r p ; xe i s refers to the number of pixels involved in 
the measurement of a particular bandpass. The results 
presented in Fig. || were obtained employing data and er¬ 
ror spectra in which Eq. (|l9]) worked properly. However, it 
is straightforward to see that, in this expression, a simul¬ 
taneous combination of large and small <r(A) values would 
lead to a poor agreement between both terms. 

For example, discrepancies between Gonzalez’s formu¬ 
lae and numerical simulations are apparent in the mea¬ 
surement of the Mg 2 index in the spectra of the central 
dominant galaxy of the cluster Abell 2255 (observed with 
the TWIN spectrograph at the 3.5m Telescope of Calar 
Alto, August 1994), as it is shown in Fig. |j. These differ¬ 
ences are due to the fact that a bright sky-line falls within 
the central bandpass of the Mg 2 index. Sky subtraction 
during the reduction process introduces a larger error in 
the pixels where sky-lines are present. As a result, the 
reduced error spectra exhibit, simultaneously, pixels with 
very different <r(A) values, and Eq. (|l9|) is a poor approx¬ 
imation. 


6. The new formulae 


We have derived a more accurate set of analytical formulae 
to compute errors in line-strength indices. The expected 
random error in the measurement of an atomic index is 
given by 


<y [4] = cr 



s( A) 

C(A) 


dA 


( 20 ) 


In practice, this integral must be transformed into the 
summation 


o' [4] - & 


N pixels 

°E 


i= 1 


S(\i) 
C( Ai) 


= 0 it [2] 


( 21 ) 


where 0 is the dispersion (in A/pixel), assuming a lin¬ 
ear wavelength scale, and IVpixeis the number of pix¬ 
els covering the central bandpass (note that, in general, 
fractions of pixels must be considered in the borders of 
the bandpasses). I is a function of 2 x Np- lxe \ s variables 
!{■■■, S(Xi ),..., C(Xj), ...) which verify 

cav(S{Xi),S(Xj)) = 0 Vi,j e [1, Arpi xe i s ], i ^ j {22) 

cov(S(Xi), C{Xj)) = 0 Vi, j g [1, Apixeis] (23) 

cov(C(Aj), C(Xj)) A 0 Vi, j G [ 1 , Apixeis], i ¥= j (24) 

since C(A) is computed from Eq. ©■ Taking this result 

into account: 



Relative Error (simulations) 


Fig. 3. Comparison of the Mg 2 relative errors employing 
Gonzalez’s formulae and numerical simulations. Filled trian¬ 
gles refer to a sample of 40 bright stars from the Lick Li¬ 
brary, whereas open triangles correspond to Mg 2 measurements 
along the radius of the cD galaxy in Abell 2255. When using 
Gonzalez’s expressions, there is a clear underestimation of the 
relative errors in the galaxy due to the subtraction of a bright 
sky-line in the Mg 2 central bandpass during the reduction pro¬ 
cess. 



-^pixels ^pixels 

+ E E 

i=i 


(- 


dl 


(■ 


dl 


\dC{Xi)) \dC(Xj) 


cov(C(A,), C(Xj)) 


N pixels 

= E 


i= 1 


C 2 {Xi) cr 2 {Xi) + S 2 {Xi) O'c^Xi) 
C\Xi) 


iW AWe. S U S{K) ^ ( s{Xj) ^ ^ ^ ^ 

§ U 2 (Ai)J C0V( ( ( j)) . 


After some manipulation the covariance terms are 
cov(C(Xi),C(Xj)) = 


(25) 


(C(Xi)C(Xj)) - (C(Xi))(C(Xj)) = (26) 


a 2 [l\ = 


Ai (Ts b 4“ A .2 cov(S&, S r ) H- A 3 cov(5 r , Sb) H - A 4 & g r (27) 
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where we have defined the following four parameters 

Ai = 

A 2 = 

A 3 = 

A 4 = 


(A r 

Ai) (A r 

Aj) 


(A r — Af ,) 2 


(A r 

— Af)(Aj - 

A b ) 


(Ar - A b ) 2 


(Ai 

Afr)(A r 

A,) 


(Ar — A b) 2 


(A* 

~ X b )(Xj - 

A b) 


(A r - Ah ) 2 

Since St, y S r are not correlated, we obtain 

cov(C(Xi), C(Xj)) = Ai a 2 Sb + A 4 a 2 Sr 

and finally 


l [h 


N, 


pixels 


0 2 


= E 


C 2 (X i )a 2 (X i ) + S 2 (X i )a 2 c{Xi 

C\x~) 


(28) 

(29) 

(30) 

(31) 

(32) 


N pixels N pixels 

+ E E 

* =1 

where 


' SjX^SjXj) 
_C 2 (Xi) C 2 (Xj) 


(Ai &s b + A4 cr| r ) 


(33) 


a 


2 

C(A 4 ) 



A r A b 


2 4 ,+ 


f A i \b 

\ A r A b 



(34) 


2 -^pixels (blue) 

= g ^ < 35 > 

^2 ^pixels(red) 

A = (y—yyr g < 36 > 

Errors in the molecular indices are calculated through 

= (37) 

which is identical to Eq. (ED¬ 
ITsing the new formulae, the computed index errors for 
the samples presented in Figs. [| and |] completely agree 
with the results from numerical simulations (see Fig. |4). 

The errors in the 4000 A break, defined in Eq. (||), can 
be computed as 


°’ 2 [D 4 ooo] — 




2^.2 


b u JF r 


^b 


(38) 


where 


-^pixels 

r P = e E [ A ? s(ao] 

2=1 


(39) 


N pixels 

°% p = 02 E l x t ^( A *)] ( 4 °) 

2=1 

(the subscript p refers indistinctly to b or r). Note that in 
this case Tb and T r are not correlated. D 4 ooo errors com¬ 
puted in this way show a perfect agreement with numerical 
simulations. 



Relative Error (simulations) 


Fig. 4. Comparison of the Mg 2 relative errors employing the 
new formulae and numerical simulations for the samples shown 
in Fig. The agreement between both methods is complete. 


7. Estimation of S/N ratios 

As it is apparent from Fig. |], there is a clear correla¬ 
tion between the measured relative error and the (S/N- 
ratio)/A. Although such a correlation arises naturally, 
quantitative estimates of errors as a function of (S/N- 
ratio)/A are also expected to depend on index values. 

We have shown that the use of approximate formu¬ 
lae for the computation of errors can lead to misleading 
error estimates. However, performing appropriate simpli¬ 
fications, it is possible to derive simple expressions to es¬ 
timate the absolute index error as a function of the mean 
(S/N-ratio)/A. It can be shown that 


cr[Q 


Ci — C 2 I& 

sn(A) 


where, to simplify, we have defined SN(A) as 


SN(A) 


1 S{Xj) 
N Ve ^ cr(A i) 


(41) 


(42) 


(the summation extends over the three bandpasses, i.e. N 
pixels). The two constants C\ and C 2 are defined as follows 


Cl = AA C C2 


(43) 
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being AA b, AA C and AA r the bandpass widths. 
Similarly 


c[/ m ] 


C3 

sn(A) 


where 


c 3 = 2.5 c 2 log 10 e 


(45) 


(46) 


Equations © and (|fi|) can be easily employed to pre¬ 
dict the required (S/N-ratio)/A to achieve a fixed index 
error. In Fig. we show the predictions of these expres¬ 
sions for some particular indices compared with the actual 
error measurements (from Eqs. [hi] and [3?]) in a star sam¬ 
ple. It is interesting to note that the absolute error of 
molecular indices does not depend, in a first approxima¬ 
tion, on the absolute index value, although the contrary is 
true for the atomic indices. This is the reason why a larger 
scatter is apparent in panel (a), where / a in Eq. (© has 
been replaced by its corresponding arithmetic mean in the 
sample. Numerical values for ci, c 2 and c 3 , which obviously 
depend on the considered index, are given in Table 1. 

Following the same procedure with the D 4000 index: 


cp^ooo] 


D 4000 / 1 1 

V20oysN(A)" sn(A)* 


(47) 


where SN(A)^ and SN(A) r , the mean (S/N-ratio)/A in 
the blue and red band respectively, will attain, in general, 
different values (given the large wavelength coverage of 
the break). In this case, the relative D 4000 error does not 
depend on the absolute D 4000 value. 


8. Summary 

We present a set of analytical formulae to compute ran¬ 
dom errors in the measurement of line-strength indices. 
The new expressions constitute an improvement compared 
with the approaches followed in previous works, providing 
a straight method to attain reliable random error esti¬ 
mates. In particular, the derived equations should be pre¬ 
ferred to methods involving either numerical simulations 
(a computer time demanding method) or multiple obser¬ 
vations (usually prohibitive at low light levels). As a use¬ 
ful tool for observation planning, we also provide simple 
recipes to estimate the required signal-to-noise ratio to 
achieve a desired index error. 

We want to stress here that full benefit from these for¬ 
mulae can only be obtained after a proper treatment of the 
error propagation throughout the data reduction. Apart 
from this, the analytical method followed to derived the 
expressions given in Section 6 is exact, in the sense that it 
does not involve any approximation (exception made for 
the unavoidable data sampling). 



20 50 100 200 500 


(S/N—ratio)/A 


Fig. 5. Absolute atomic (panel a) and molecular (panel b) 
line-strength errors measured in the sample of 40 stars from 
the Lick library as a function of the mean (S/N-ratio)/A. The 
predictions of Eqs. © and © are plotted as dashed lines. 
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